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ABSTRACT 



A fundamental Doppler-like but asymmetric wave effect 



^ that shifts received signals in frequency in proportion to 
Q their respective source distances, was recently described 
as means for a whole new generation of communication 
technology using angle and distance, potentially replacing 
TDM, FDM or CDMA, for multiplexing. It is equivalent to 
' (-H 'wave packet compression by scaling of time at the receiver, 
O fonvertins path-dependent phase into distance-dependent 
^ shifts, and can multiply the capacity of physical channels. 

The effect was hitherto unsuspected in physics, appears to 
^^be responsible for both the cosmological acceleration and 
^ \he Pioneer 10/11 anomaly, and is exhibited in audio data. 
' ^ ^^^^ paper discusses how it may be exploited for instant, 
"^^.^assive ranging of signal sources, for verification, rescue 
g.nd navigation; incoherent aperture synthesis for smaller, 
I ^H vet more accurate radars; universal immunity to jamming 
or interference; and precision frequency scaling of radiant 
energy in general. 



> 



I. Introduction 



^ A previously unsuspected result of the wave equation, 
^ enabling a receiver to get frequency shifts in an incoming 
^ signal in proportion to the physical distance of its source, 
was recently described [1]. Its main premise, that signals 
O from real sources must have nonzero spectral spread, seems 
j> fo be well supported by the cosmological and the Pioneer 
10/11 anomalous acceleration datasets, and a remarkably 
*^ large gamut of terrestrial mysteries that can be mundanely 
resolved. This perceived support actually concerns a further 
inference that the mechanism could occur naturally in our 
instruments on the order of magnitude of 10~'** s~^ This is 
a decay corresponding, as half-life, to the age of the solar 
system, and is too small for purely terrestrial applications. 
The wave effect is consistently demonstrated with acoustic 
samples, however, testifying to its fundamental and generic 
nature. While a general exposition deserves to be made in 
due course in a physics forum, it presents fundamental new 
opportunities for intelligence and military technologies. 

Section |ll] contains a brief review of the theory of the 
effect, showing that all real signals necessarily carry source 
distance information in the spectral distribution of phase. 



analogous to the ordinary spatial curvature of wavefronts, 
and that by scanning this phase spectrum, each frequency 
in a received waveform comprising multiple signals would 
be shifted in proportion to its own source distance and the 
scanning rate. Section |lll]presents the general principles of 
realization by both spectrometry and digital means. Section 
ITVl shows how this enables separation of signals by physics 
instead of modulation, whose implications for information 
theory have been discussed [ibid.]. Intelligence and military 
possibilities are considered for the first time in Section jV] 

II. Rubble's law shifts, temporal parallax 

The sole premise for the effect, as mentioned, is the 
nonzero bandwidth of a real signal. The Green's function 
for the general wave equation concerns an impulse function 
5{x,t) as the elemental source. Its Fourier transform is 



F[5(0] = / 5{t)e-''"dt = l 



(1) 



which says that all spectral components of an impulse start 
with the same phase. The source itself is the only common 
reference across any continuous set of frequencies, hence 
the spectral phase contours indicate its distance (Fig. [l]). 




distance from source 



Fig. L Phase contours and gradient for a source impulse 

The main difficulty in tapping this available source of 
source distance information is that it lies in differences of 
phase across incoming frequencies, and phases are hard to 
measure accurately. However, if we scanned the spectrum 
at rate dk/dt, then, for each contour (|) and source distance 
r, we should encounter an increasing or decreasing phase 
proportional to r and the integration interval Ak, as shown 
by the shaded areas in the figure. The trick that yields this 
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remarkable asymmetric effect is in the definition that a rate 
of change of phase is a frequency, so that the instantaneous 
measures of the spectral phase contours obtained from the 
scanning have the form of frequency shifts proportional to 
the slopes, and therefore to the respective source distances. 
These distance-dependent shifts resemble Hubble 's law in 
astronomy but are strictly linear due to their nonrelativistic, 
mundane origin, and are rigorously predicted by basic wave 
theory, as follows. The instantaneous phase of a sinusoidal 
wave at (r,f), from a source at r = 0, would be 

(|)(k,(o,f) = k-r-(Of , (2) 

and leads to the differential relation 

A(|)|^_, = A(cL)f)+k-Ar + Ak-r . (3) 

The first term A((0 t) clearly concerns the signal content 
if any, and has no immediate relevance. The second term 
k • Ar expresses the path phase differences at any individual 
frequency, and is involved in the ordinary Doppler effect, as 
it describes phase change due to changing source distance, 
as well as image reconstruction in holography and synthetic 
aperture radar (SAR), as the reconstructed image concerns 
information of incremental distances of the spatial features 
of the image. The third term, clearly orthogonal to other 
two, represents phase differences across frequencies for any 
given (fixed) source distance. Then, by varying a frequency 
selection k at the receiver, we should see an incremental 
frequency, or shift. 



distance 
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(4) 



(Note that k can only refer to the incoming wave vector 
in the denominator.) As this holds for each individual O), 
from equation ([2]), the signal spectrum would be uniformly 
shifted by the normalized shift factor 



5(0 _ 1 3(|) 
(0 m dk 



with 



(5) 



dk r dk 
' dt kc dt 

^ = k~\dk/dt) and a = P/c , 

where t denotes time kept by the receiver's clock. Further, 
the shifts reflect only the instantaneous value of P, which 
is the normalized rate of change at the receiver. 

Fig. |2] shows that 5(0 is a temporal analogue of spatial 
parallax, with the instantaneous value of a = P/c serving 
in the role of lateral displacement of the observer's eyes: 
at each value tti or a2, the signal spectrum f (co) shifts to 
^((Oi) or !F ((O2), respectively, and if the distance to the 
source were increased to r+5r, the spectrum would further 
shift to !F ((O3). Unlike with ordinary (spatial) parallax, the 
receiver can be fixed and monostatic, and yet exploit the 
physical information of source distance. 




Fig. 2. Temporal parallax 



III. Physics of realization 

There are three basic ways of accomplishing Fourier 
decomposition or frequency selection in a receiver: using 
diffraction or refraction, a resonant circuit, or by sampling 
and digital signal processing (DSP). Realization of the 
shifts in all three methods has been described, along with 
how this "virtual Hubble" effect had remained unnoticed so 
long [ibid.]. A review of the diffraction and DSP methods is 
necessary as theoretical foundation for the application ideas 
to be discussed. These methods would also be applicable 
for array antennas and software-defined radio, respectively. 

A. Diffractive implementation 

Diffractive selection of a wavelength k concerns using a 
diffraction grating to deflect normally incident rays to an 
angle corresponding to the selected k, as determined by 
the grating equation nX = ZsinS, where n denotes the order 
of diffraction, and / is the grating interval. 

The property exploited is that rays arriving at one end of 
the grating combine with rays that arrived at the other end 
a little earlier. If we could change the grating intervals I 
inbetween, the rays that get summed at a diffraction angle 
would correspond to changing grating intervals l{t), as 
depicted in Fig. [3] Their wavelengths must then relate to the 
grating interval as nd'kjdt = (dl/dt) sin0. Upon dividing 
this by the grating equation, we obtain the modified relation 
X-\di/dt)=T-\d7/dt) = -^. 



(rays at ^o) 
grating 



(rays at t^) 1 lens 




Fig. 3. Time- varying diffraction method 
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This result is not merely a coincidence of different 
wavelengths at the focal point, as would be obtained with 
a nonuniform grating. Each of the summed contributions 
itself has a time-varying X, so the variation over the sum is 
consistent with the component variations. A grating gathers 
more total power than, say, a two-slit device, and the waves 
arrive with instantaneously varying phases 

^ = -ro = ^ + Vr^-r + Vr4-r' + Vk^-k , (6) 
at at 

where the first term d(^/dt is the signal (cof) contribution 
in equation Q and is equal to -co; the second term (in r) 
concerns the real motion and Doppler effect if any, to be 
ignored hereon; the third term (in r') represents a similar 
Doppler shift from longitudinal motion of either end of the 
receiver, which would be generally negligible for r' <^ r; 
and only the last term concerns equation 

B. DSP implementation 

In DSP, k's, are determined by the sampling interval T , 
and can therefore be continuously varied by changing T . 
The discrete Fourier transform (DFT) is defined as 



N-\ 



«=o 



(7) 



with the inverse f{nT) 



im(£iTT 
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F{m(>ij) , 



m=0 



where / is the input signal; T is the sampling interval; 
is the number of samples per block; and 0)7- = 2% /NT. The 
inversion is governed by the orthogonality condition 



N-l 



n=0 



^im(xijT iikr—ntHjT) 



1 



\ — gi{m-n)/N 



(8) 



where 5„,„ = 1 if m = n, else 0. These definitions as such 
suggest that the instantaneous selections (Or = kc can be 
varied via the sampling interval T. 

To verify that a controlled variation of T will indeed 
yield the desired shifts, observe that in equation Q, both 
the real Doppler terms, in r and r', can be ignored, as r' <C r 
for sources of practical interest. The surviving term on the 
right is (d(^/dk){dk/dt) where d(^/dt = —0), as before, and 

die _ 1 dd)T _ 1 271 \ _ 271 dT _ -^^ dT 
'dt ~ 'c'dF ~ 'cdt [nT J ~ ~ NcT^~d7 ~ ~ T~dt ' 

so that, corresponding to equation Q, we do get 

k-\dk/dt) = p = -T-\dT/dt) , (9) 

confirming the desired effective variation of k's. □ 

Fig. |4] explains the result. The incoming wave presents 
increasing phase differences 5(|)i, 5(|)2, 5(|)3, ... within the 
successive samples obtained from the diminishing intervals 



5Ti = Ti — To, 572 = T2 — Ti, 5T^ = Tt, — T2, etc. From the 
relation (Oj = kc = 2k /NT, gradients of the spectral phase 
contours can be quantified as 

dk 271 
df~Mdf~~'N^M 
so that by equation each phase gradient reduces to 

3(|) dT -271 1 dT (Or 1 dk dk 



dT dt NcT dk T dt c dk k dt dk dt 
identically, validating the approach. 



(10) 



incoming wave 




To sampling Ti T2 T3 T4 T5 

Fig. 4. Variable sampling method 

C. General operating principles 
For a desired shift z at a range r, equation (|5]l yields 



hk 
~k 



5r 
r 



cz 5t 



(11) 



for normalized incremental change of the receiver's grating 
or sampling interval required over a small interval dt. In a 
DSP implementation, appropriate at suboptical frequencies, 
this determines the needed rate of change over a sampling 
interval 5T. As 5T would be nominally chosen based on the 
(carrier) frequencies of interest and not on the range, the 
nominal rate of change for achieving a useful z at a given 
range r would be independent of the operating frequencies. 

For example, for a 1 GHz signal, a suitable choice of 
sampling interval is 5t = 100 ps, regardless of the range. 
Chosing z = 2, which is fairly large in astronomical terms 
but convenient for the present purposes, we obtain dk/k = 
6 X 10~^/r per sample (taking c = 3 x 10** m/s). At r = 1 m, 
the result is a somewhat demanding 6 x 10~^, or 6%, per 
sample, but for r = 1, 10 or 1000 km, it is 6 x 10"^, 6 x 
10~^ and 6 x 10~'' per sample, respectively. 

To be noticed is that since equation ([TT]) prescribes a 
normalized rate of change, the instantaneous rate of change 
must grow exponentially in the course of the observation, 
since by integrating equation ([TT]) . we get 



^k = AT- 



gCzAt/r 



(12) 



where At denotes the total period of observation. While an 
exponential variation is generally difficult to achieve other 
than with DSP, two other problems must be also contended 
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with: the total variation over an arbitrarily long observation 
would be impossible anyway and it could exceed the source 
bandwidth even over a fairly short observation. 

The only solution is to split the observation time into 
windows, as in DFT, resetting the sampling interval and 
its variation at the beginning of each window, so that the 
variation is only continuous within each window. 

For example, the normalized rate 6 x 10"*^ calculated 
above for 100 ps sampling and 10 km range amounts to 
3.773 X 10^^" for each second of observation! Over a 1 (j.s 
window, however, the sampling rate would change only by 
the factor 1 .0006. This is also the spectral spread required 
of the wavepackets emitted by the source, and conversely, 
the window can be set to match the spread. 

It is possible to cascade multiple stages of gratings or 
DSP, in order to multiply the magnitude of the shifts. 
Cascading seems straightforward with optical gratings, but 
in DSP, as the data stream is already broken up into discrete 
samples by the first stage, and the sampling interval is to be 
varied in each successive stage relative to its predecessor, 
it becomes necessary to interpolate the samples between 
stages. Interpolation is similarly needed for the reverse shift 
stage in distance-based signal separation (Section IIV-BI ). 

These shifts have been verified only by simulation for 
electromagnetic waves. For sound data (c f» 330 m/s), the 
software does reveal nonuniform spreading of subbands of 
acoustic samples, consistent with the likely distributions 
of their sources. Testing in an adequately equipped sound 
laboratory, as well as with radio waves, is still needed. 

D. Differential techniques 

The effect thus scales well from terrestrial to planetary 
distances, but the incremental changes required per sample 
are very small. Errors in these changes, from inadvertent 
causes including thermal or mechanical stresses, could lead 
to large errors in distance measurements using the effect. 
The linearity of the effect can be exploited to offset such 
errors by differential methods. From equation dD, we have 

Az = aAr + rAa + o(Ar,Aa), (13) 

relating first order uncertainties. First order error due to an 
uncertainty Aa can be eliminated by measurements using 
different a and using the differences. 

Specifically, an error Az would limit the capability for 
source separation, presuming r is known. By designing a 
receiver to depend on the difference between two sets of 
shifts for the same incoming waves, e.g. by applying the 
DSP method twice with different values of a, a Az error can 
be eliminated. Conversely, for ranging applications, where 
z is the measured variable, from r = z/a, we get 

Ar = Az/a-rAa + o(Az,Aa), (14) 



so that the first order error Ar can be once again eliminated 
using differences. This is a temporal form of triangulation, 
illustrated by the lines for tti and a2 that converge at r in 
Fig. |2] Higher order differences would yield more accuracy. 

E. Phase shifting and path length variation 

The main difficulty with the method of Section IIII-AI is 
ensuring uniformity of the grating intervals even while they 
are changing. The method is not realizable by acousto-optic 
cells for this reason. Controlled nonuniform sampling and 
sample interpolation pose similar difficulties for DSIo. 

The solution is to instead vary the path length after a 
fixed grating, say using the longitudinal Faraday effect and 
circular polarization. As the variation now concerns a bulk 
property, it would be easier. The proof is straightforward. 
The analogous simplification for DSP is to modify, instead 
of the sampling interval T, the forward transform kernel to 

yv-i 

f(")^«o)= 1^^^™"'''^/K), (ait) = (aoe^"^, (15) 

n=0 

i.e. apply changing phase shifts to the successive samples 
corresponding to the path length. This is more complex, but 
avoids interpolation noise and access to the RF frontend. 

IV. Source separation theory 

A. Transformation kernel 

Fourier theory more generally involves the continuous 
form of the orthogonality condition ([8]l, given by 

where 8() is the Dirac delta function. All of the methods 
described require continuously varying the mechanism of 
frequency selection. This is equivalent, in a basic sense, to 
varying the receiver's notion of the scale of time, the result 
being a change of the orthogonality condition to 

Jt Jr (17) 

= e''"'^5[(i>A-0)] , 

where A = A(r) = (1 +OCr) from equation ([5]), but is also 
equivalent to A{t) = (1 +act), via the relation c = r/t. 

How did we get to equation (fTTl)? Notice that equation 
([T6l ) already contains the frequency selection factor e"^' - 
the A in the exponent is the variation of (O provided by the 
methods of Section |llll The second factor must 
correspondingly represent the incoming component picked 
by the modified selector; its only difference from equation 
([T6l ) is the A multiplying the path contribution to phase, kr. 

'Radio telescopes, for instance, incorporate 1 or 3 bit sampling at an 
intermediate frequency for correlation spectroscopy. Interpolation would 
add more noise to this already low phase information. 
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Why should this A multiply kr and not (Of? The answer lies 
in the basic premise of equations ©-([S]! that our receiver 
manipulates the phase of each incoming Fourier component 
individually, hence COf, representing the signal component 
of the instantaneous phase, is not touched. 

The kr term is affected, however, since the summing 
builds up amplitude at only that value of k which matches 
65 = kc. This aspect is readily verified by simulation, and it 
defines a lower bound on the window size (Section IIII-CI ) 
in the optical methods, since the window reset would break 
the photon integration in a photodetector. 

Equation ([TT] ) is more general than equation (fT6l ). as the 
latter corresponds to the special case of A = 1 or a = 0. 
This property has a fundamental consequence impacting all 
of physics and engineering: Since a refers to the physics of 
the instruments and could be arbitrary small but nonzero, 
there is fundamentally no way to rule out a nonzero (X in 
any finite local set of instruments. The only way to detect 
a nonzero a is to look for spectral shifts of very distant 
objects as a large enough r would yield a measurable z{r^ 

B. Distance division multiplexing™ 

Almost all the applications to be discussed critically 
depend on the fundamental separation of sources enabled 
by the present wave effect. Using the notation of quantum 
theory, we may denote incoming signals by "kets" |), and 
the receiver states as "bras" (|, and rewrite equation (fT6l ) 
as ((o|a),r) = e'*^''((o|a)). Equation (fTTl) then becomes 



CO 



JkrA(r) 



5 ( CO 



CO 



(18) 



where (co,Jco/Jf| and (co| are the modified and original 
states of the receiver, respectively. The sampling clock and 
grating modifications then correspond to the operator H 

CO , 



H\(0,r) 



'A(r) 



(19) 



^ A personal hunch that something hke this could be responsible 
for the cosmological redshifts, had prompted informal prediction of the 
cosmological acceleration to some IBM colleagues in ca. 1995-1996. 

The modified eigenfunctions are known in astrophysics as 

the photon and particle eigenfunctions over relativistic space-time [2]. 



The natural value of [a ~] 10 



cited in the Introduction can come 



from the probability factor e ^ I^bT cumulative lattice dislocations 
under centrifugal or tidal stresses causing creep. At T ~ 300 K, it is 
3 X 10"" s"' at IV = 1 eV, 1.5 x 10""* s"' at IV = 1.7 eV, etc. The 
point is that solid state theory mandates such creep, but it's untreated 
in any branch of science or engineering. Many of these details were put 
together in mundane arxiv.org articles [4], to be eventually compiled 
into a comprehensive paper. Thanks to NASA/JPL's continued portrayal 
of the anomaly as acceleration, all other explanations tested or offered 
have been more obvious or exotic, and totally fruitless (cf. [3]). 



of an incoming wave state |co,r). These equations attribute 
the shift to the incoming value CO, instead of to CO selected 
instantaneously since for the operator formalism, (co| must 
represent an eigenstate resulting from observation. 

Equations dSjl and ( [T9l ) both say that as in the Doppler 
case, the shift is proportional at each individual frequency 
CO. This also means that the spectrum expands by the factor 
A, as illustrated in Fig. |5] for the case of a common signal 
spectrum f (co) emitted by two sources at distances ri and 
^2 > ri , respectively. For a rate of change factor a applied 
at the receiver, the signals will shift and expand by Ai = 
(1 +ari) and A2 = (1 + ocr2) > Ai, respectively. 



A2r2 




Gi 



frequency 
Fig. 5. Separation by source distance 

Despite the increased spreading of the spectrum, there 
is opportunity for isolating the signal of a desired source if 
its shifted spectrum comes out substantially separated from 
its neighbours, as shown. We could, for instance, apply a 
band-pass filter Gi to the overall shifted spectrum HF\ + 
HF2, such that Gi{HFi+HF2) ^ GyHFy. Writing // as a 
function of a, by equations (fTTl ) and (fT9l ). we get 



//-1(a) = -//(a) 



(20) 



We could apply a combination of mixing and frequency- 
modulation operations to shift and compress GiHFi back 
to f«/^i(co), or just apply a second H with a negative a. In 
general, a set of distance-selecting "projection operators" 
H~^GiH, can thus be defined by the conditions 

H-^GiH ^F(r„co) PS F(r/,co) or H'^GiH ^ 5,y (21) 



With base-band prefilters G/i^ ^ Fi, they yield 
H-^GM = hijGj so that GiH = HGi 



(22) 



As H is parametrized by a independently of /, HGiH 
provides spectral separation without prior knowledge of r,. 
It is also independent of signal content, wherein the source 
coordinates could be supplied or included via modulation, 
but since H depends only on the path phase contribution, it 
provides separation at a more basic level. Coordinate-based 
spread spectrum coding could be combined, for example, to 
differentiate sources at the same distance, as an alternative 
to or to improve over current phased array antennas. 

It has been further suggested [1] that DDM raises the 
theoretical capacity of a channel to ~ 2CL/X » C, where 
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C is the traditional Shannon capacity for a source-receiver 
pair using the channel, by admitting additional sources into 
the channel at intermediate distances. Further, if combined 
with directional selectivity say from phased array antennae, 
we would get true space-division multiplexing {ibid.^. 

V. Military and intelligence applications 

The above theory represents a long overdue rigorous 
analysis of the spectral decomposition of waveforms by a 
real, and therefore imperfect, instrument. The nature of the 
error was inferred, as mentioned, from an immense gamut 
of astronomical, planetary and terrestrial data, spanning all 
scales of range ever measured. A natural occurrence of the 
effect subsequently deduced from creep theory (footnote, 
page [5ll could have been contradicted in at least one set 
of data, but remains totally consisten{^. Instead, we clearly 
face a converse problem in physics itself, that its patterns 
have been so well explored on earth that only astronomical 
tests can expose remaining shortcomings [5]. The known 
laws of physics only reflect understood mechanisms, like 
interference and the Doppler effect and only partiall}]^. We 
once again face a need to lead physics by engineering 

A. Precision, power-independent frequency shifting 

This application would itself be a direct, visible test of 
the effect, the idea being that the method could be applied 
also to a proximal source at a precisely known distance in 
order to uniformly scale its spectrum with great accuracy. 
Accuracy is expected because the source distance can be 
accurately set by simple mechanical means, independently 
of the exponential control of a. The mechanism would be 
the first means for shifting frequencies that is independent 
of the nature and energy of the waves, and their frequencies 
and bandwidth, yet realizable by a mundane, static means. 
The main difficulty is the magnitude of a needed for source 
distances of under 1 m, but it should be eventually solvable. 

Likely applications include transformation of high power 
modulated GHz carriers to THz or optical bands, and 

■'Consistency with GPS datasets has also been recently verified. 
Calculations from GPS base stations data indicate an ongoing rising of 
land at a median rate of 1.6 mm/y (see http://ray .tomes .biz), 
large enough to easily accommodate a nontectonic apparent expansion 
of the earth at Hq x 6.371 km ~ 0.437 mm/y, for the natural occurrence. 

^E.g., diffractive corrections in astrophysics and quantum field the- 
ories are limited to the Fraunhofer and Fresnel approximations, both 
being limited, by definition, to total deflections of < 7l/2. Cumulative 
diffractive scattering would contribute, as known in microwave theory, a 
decay e^"'' to the propagation. This makes all wavefunctions in the real 
universe Klein-Gordon eigenfunctions, and gives light a rest-mass... 

^Recalling the classic case of thermodynamics. The present effect 
implies that photons, traditionally viewed as indivisible and immutable, 
are reconstituted by every real, imperfect, receiver. Though it is yet to be 
demonstrated with light, the theory of Sections |ll] and |lll] is fundamental 
and does not permit a different result for quanta. 



efficient, tunable visible light, UV or even X- or y-rays, all 
with controlled coherence and without nonlinear media. 

B. Monostotic ranging and passive radar 

This application was immediately envisaged when the 
effect was first actually suspectecj^, taking from the well 
known notion of the cosmological distance scale available 
from the Hubble redshifts. By providing a "virtual Hubble 
flow" view that can be activated and scaled at the receiver's 
discretion, the present methods enable ranging, or distance 
measurement, of any source that can be seen or received, at 
only half the round-trip delay, and zero transmitted power. 

To compare, as traditional active radar depends on round 
trip times (RTT), its range is fundamentally limited by 
the transmitter power. While Venus and Mars have been 
explored by radar, for example, the ranging of other planets 
and of astronomical objects in general is largely limited 
to spatial triangulation, including using the earth's orbit 
itself as the baseline. Tracking of orbiting satellites places 
similar demands on the radar transmitter power P, as the 
range R is governed by the "radar power law" P oc R^. This 
makes current (active) radars generally bulky. The power 
law is relaxed for cooperative targets with transponders, to 
P<x R^, where P denotes transponder power, and the range 
is obtained from the RTT of a transponded signaj^ The 
present methods make this lower power law available for 
all targets, and also eliminate the need for a spatial baseline 
for parallax, as the baseline is now given by the range of 
variation of a at the receiver. As shown in Section ITlI-CI a 
very large baseline seems to be indeed possible with DSP 
alone for both terrestrial and earth orbit distances. 

Current passive radar technology, like Lockheed's Silent 
Sentry, involves coherent processing of the received scat- 
ter of energy from television broadcast (see US Patent 
3,812,493, 1974) and cellular base stations. The accuracy 
is again dependent on information representing overall trip 
times, and hence entails coherent processing. 

Accurate ranging is now possible via triangulation using 
temporal parallax, as explained in Section IIII-DI The 
difference is that coherent processing can provide imaging 
with subwavelength precision, but triangulation would be 
simpler for ranging and tracking. We still need to illuminate 
silent targets, but the processing could be simplified for the 
same accuracy. The accuracy of tracking could improve for 
radiating targets given the lower overall trip time. We can 
also use temporal triangulation to simplify traditional radar 
and as a fast, efficient cross-check. 

^The idea of such an effect had been disclosed, by appointment, to an 
IP attorney on the morning of 2001.9.11. The mechanism itself (Section 
Hit was uncovered only much later in 2004. 

^This is generally what's used for deep space probes [7], [8]. 
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C. Overcoming interference, jamming and noise 

The capabiUty for separating sources regardless of the 
signal content or modulation is also a guarantee that a 
desired signal can be isolated from any interfering signals 
at the same frequencies and bearing the same modulations 
or spread-spectrum codes, so long as the interfering sources 
are at other directions or distances from the receiver. 

Simulation shows that even signals in destructive inter- 
ference can be separated. The criterion for separation is 
simply the bandwidth to distance ratio: denoting the low 
and high frequency bounds of the signal bandwidth W by X 
and 9{ , respectively, each of equations (|5]l or (fTTl ) implies 



{\+ari)9{ < (1 +ar,+i)i: 



(23) 



as the condition for separation between the ith and (/+ l)th 
sources, assuming r, < r,+i (see Fig. [5]l. This means 



W 



{5ri/ri)L-W 



5r, L 
7~W 



1 



(24) 



where W = ^ — L and 5r, = r,+i — r,. As separation is 
only possible for a > 0, we need to have 5r, > r, W /l. It 
then follows, however, that sufficiently narrow subbands of 
the total received signal would be "source-separable" even 
when the condition should fail for W as a whole. 

All of the limitations are thus purely technological, e.g. 
how many and how narrow subbands we can construct, 
how linear we can make the subband and band-pass (G) 
filters (since phase envelope distortions will alter the shift), 
and how large the stop-band rejection we can get. The stop- 
band rejection ratio is important in overcoming jamming, 
and the separation as such should suffice, in principle, for 
filtering out truly extraneous noise like lightnings. 

D. Incoherent aperture synthesis 

In a conventional SAR, multiple targets or features are 
imaged from the echos received on a moving platform, 
typically an aircraft or satellite, as depicted in Fig. [6] 
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Fig. 6. Aperture synthesis 



If the onboard oscillator could maintain coherence for 
the duration of the flight, i.e. not suffer random phase 
shifts, an entire target region could be imaged by a simple 
Fourier transform, after correcting for the Doppler shifts in 
the received echos due to the radar's own motion. At each 
point Q on the flight path, echos are received from multiple 
features A, B, etc. and corresponding to multiple previous 



locations Pi, Pi, ... of the transmitter. Fourier techniques 
and their optical holographic counterparts are most efficient 
for unravelling the information. 

As stated in Section IV-BI the fundamental reason for 
coherent processing is our current dependence on RTT of 
radar echos for the range information. The dependence is 
obviated by the wave effect. In the example scenario, echos 
received at Q from target features differing in direction, like 
A and B, can be separated by the phase information from 
an antenna array. The echos from the same direction, such 
as from A and D, can now be similarly separated using the 
methods of Section |lll] - in particular, the changing delay 
transform of equation (fTSl) can be applied to recordings of 
the received echo waveform. 

Coherent processing ordinarily has two advantages that 
would seem to be impossible in any other approach: it 
tends to curtail noise and is generally precise to within a 
wavelength. We would get suppression of external noise, as 
discussed in Section IV-Cl and subwavelength resolution, as 
shown in Section ITlI-Di so that applications including hand- 
held short-range carrier-deck operation, for which coherent 
processing is difficult, become feasible. 
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